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Glossary 

Acceleration mechanism: Any process, or sequence of processes, that results in a small fraction of the particles in 
a plasma each gaining an energy greatly in excess of the thermal energy. 

Heliosphere: The region around the Sun populated by plasma escaping from the Sun in the solar wind, and extending 
to well beyond the orbits of the planets. 

Magnetosphere: The region around a planet or star populated by plasma on magnetic field lines tied to the planet 
or star. 

Resonance: The gyroresonance condition (6) corresponds to the wave frequency being an integral, s, multiple of the 
particle gyrofrequency in the rest frame of the particle gyrocenter; s = is also called the Cerenkov condition. 
Shock: A fast magnetoacoustic shock is a discontinuity in plasma fiow that propagates into the upstream plasma 
faster than the Alfven speed, compressing both the plasma and the magnetic field. 
Solar corona: The hot (> 10^ K) outer atmosphere of the Sun. 

Solar flare: An explosion in the solar corona that lead to copious acceleration of fast particles. 

I. BACKGROUND AND CONTEXT OF THE SUBJECT 

A typical astrophysical or space plasma consists of thermal plasma components, a magnetic field, various nonthermal 
distributions of fast particles and various kinds of turbulence involving different kinds of waves. There is no single 
mechanism by which the fast particles are accelerated, and it is now recognized that the data require a rich variety of 
different mechanisms operating under different conditions. Three generic mechanisms are emphasized in this review: 
stochastic acceleration, diffusive shock acceleration and acceleration by parallel electric fields. Before discussing these, 
it is appropriate to summarize the observational discoveries that drove the formulation of the theoretical ideas. 

The major observational discoveries in this field were almost all unexpected and initially quite baffling. The first 
notable discovery was that of cosmic rays (CRs), by Victor Hess in 1912. Hess measured the rate of ion production 
inside a hermetically sealed container in a balloon fiight to an altitude of 5300 meters, and concluded that "The results 
of my observation are best explained by the assumption that a radiation of very great penetrating power enters our 
atmosphere from above." By the late 1920s CRs were identified as charged particles, and not photons, and due to 
their isotropy it was inferred that they come from the Galaxy. Later in the 1930s it was accepted that a component 
of the CRs originates from the Sun. The next major discovery followed the development of radio astronomy in the 
1940s. This led to the identification of a wide class of sources radiating by a nonthermal mechanism, identified in 
the late 1940s as synchrotron radiation. The presence of highly relativistic electrons in synchrotron sources raises 
the question of how they are accelerated. Subsequent development of X-ray and gamma-ray astronomy placed severe 
constraints on the acceleration of the relativistic electrons in synchrotron sources. 

Space physics began in the late 1950s with the discovery of the Earth's Van Allen radiation belts. Further exploration 
led to the identification of energetic particle distributions in essentially all space plasmas. Observations from spacecraft 
provide in situ data on particles and fields, in the Earth's magnet osphere, in other planetary magnetospheres, and 
throughout the heliosphere. In the case of the magnetosphere. Earth-based data on precipitation of particles, notably 
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in the auroral zones, and on radio waves originating from the magnetosphere complement the spacecraft data on 
particles and waves in situ. 

The currently accepted idea on the origin of Galactic CRs (GCRs) is associated with supernova explosions was 
originally proposed by Baade and Zwicky in 1934 [BZ34]. Supernovae are the only adequate source of the energy 
needed to replace cosmic rays, whose lifetime in the Galaxy is about 10 ''yr [GS64]. Our present-day understanding 
of the acceleration mechanisms involved had its forerunners in two papers by Fermi [F49, F54]. In the first of these 
[F49], Fermi suggested that the acceleration of GCRs is due to reflections off moving interstellar clouds. Particles 
gain energy in head-on collisions with clouds, and lose energy in overtaking collisions; a net average gain results from 
head-on collisions being more probable than overtaking collisions. This model is recognized as the flrst example of 
what is called stochastic acceleration. In the second paper [F54], Fermi proposed acceleration between 'closing jaws' 
in which all reflections are head-on and hence energy-increasing. This idea is an essential ingredient in diffusive shock 
acceleration (DSA), proposed independently by various authors in 1977-8 [A77, K77, B78, BO]. DSA is now widely 
accepted as the acceleration mechanism for GCRs and for relativistic electrons in most synchrotron sources. 

Both stochastic acceleration and DSA require that the particles be scattered efficiently in order to be accelerated. 
Very efficient scattering is also required to explain precipitation of particles trapped in the magnetosphere. The 
problems associated with the scattering and acceleration of fast particles in astrophysical plasmas essentially defined 
a new field of 'plasma astrophysics' that emerged in the mid 1960s. The efficient scattering is attributed to resonant 
wave-particle interactions. 

A third generic mechanism is acceleration by an electric field, £^11 , parallel to the magnetic field. There is over- 
whelming evidence that such acceleration occurs, most notably in the Earth's auroral zones. The acceleration is 
associated with a parallel current; historically, the existence of such parallel currents in the auroral zones is perhaps 
the oldest idea relevant to the theory of acceleration, having been predicted by Kristian Birkeland in 1903. However, 
understanding such acceleration is problematical. Astrophysical plasmas are typically described using ideal magneto- 
hydrodynamics (MHD), which requires = 0, and theories involving are necessarily outside the scope of MHD 
MHD. Moreover, and E\\ accelerates particles of oppositely signed charge in opposite directions, setting up a current 
fiow, that is opposed by inductive effects, and the resulting charge separation tends to neutralize the E\\. Hence, 
despite its apparent simplicity, acceleration by E\\ is the least understood of the three generic mechanisms. 

II. STOCHASTIC ACCELERATION 

Second-order Fermi acceleration [F49] is now identified as an archetypical form of stochastic acceleration, and the 
term Fermi acceleration is sometimes used in a generic sense to describe stochastic acceleration. 

A. Historical remarks 

Fermi-type acceleration occurs when particles are refiected many times from a distribution of moving magnetic 
inhomogeneities. In the original version of the mechanism Fermi [F49] was concerned with GCRs bouncing off 
magnetized interstellar clouds. Such acceleration is more efficient for more frequent, smaller energy changes than for 
less frequent, larger energy changes, and this led to the recognition that a more efficient Fermi- type acceleration can 
result from MHD turbulence [T55, K56, D56, P57, PT58]. One form of the idea is that the slow compression and 
rarefaction of B, associated with the MHD turbulence, conserves the magnetic moment, p\/B^ and when combined 
with a scattering process that maintains the isotropy of the particles, this implies a 'betatron' acceleration, or 'magnetic 
pumping' [S57, B58]. A closer analogy with Fermi acceleration is acceleration due to the small fraction of the 
particles that reflect from the compressions in which propagate at the Alfven speed, vx- Besides the betatron- type 
acceleration and the reflection off moving compressions, there is also a transit acceleration [S65] due to particles 
diffusing through the compressions and rarefactions. All these energy changes contribute to stochastic acceleration 
by MHD turbulence [KF71]. 

The mathematical description of Fermi-type acceleration is in terms of isotropic diffusion in momentum space. In 
the early literature stochastic acceleration was treated using a Fokker-Planck approach that includes a term that 
describes a systematic acceleration and a term that describes a diffusion in energy. In was shown by Tverskoi [T67] 
that for Fermi's original mechanism this equation reduces to an isotropic diffusion in momentum space. The same 
equation is derived for acceleration by isotropic turbulence [S66, HS67, KF71]. However, MHD turbulence is never 
isotropic: it consists of a spectrum of MHD waves, in both the (non-compressive) Alfven mode and the (compressive) 
fast magnetoacoustic mode. 
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B. Diffusion in momentum space 

The isotropic diffusive equation in momentum space that describes stochastic acceleration is [T67, T68, KF71] 



d{f){p) _ 1 d 
dt dp 



(/>(P), (1) 



DpM = ^a'^(^-%) , (2) 



where {f){p) is the particle distribution function, /(p), averaged over pitch angle. The diffusion coefficient is 

4:V 

where ua is the acceleration rate. The final factor in (2) does not appear in simpler treatments of Fermi-type 
acceleration; this factor may be ignored only in the limit of particle speeds much greater than the Alfven speed, i.e., 
for V VA- The meaning of ua is most easily understood by estimating the mean rate of acceleration [T66, M68] 

where e = (ini?c^ +p|c^ -\-p\c^Y^'^ is the energy. In the limit v ^ vx (3) gives 

) ~ ^Apc, (4) 



dt I 

which reduces for highly relativistic particles to the familiar form {de/dt) ^ va£- The acceleration rate may be 
estimated as 



where 5B is the magnetic amplitude in the waves and uj is their mean frequency. A remarkable feature of (5) is 
that, although effective scattering is an essential ingredient in the theory, the acceleration rate is independent of the 
details of the scattering. The parameters that appear in (5) refer only to the MHD turbulence that is causing the 
acceleration. 

C. Treatment in terms of resonant interactions 

The treatment of Fermi- type acceleration by Achterberg [A81] shows that it results from a resonant damping of 
magnetoacoustic waves in the presence of efficient pitch-angle scattering. The general gyroresonance condition is 

UJ — sQ. — k\\V\\ - 0, (6) 

where (jJ is the wave frequency, k\\ is the component of its wave vector parallel to B, s = 0, ±1. ... is the harmonic 
number, O = \q\B/^m is the relativistic gyrofrequency of the particle with charge mass m, Lorentz factor 7 and 
velocity v\\ parallel to B. The effect of gyroresonant interactions on the particles and the waves is described by the 
quasilinear equations, which are written down in the Appendix. The particular resonant interaction of interest here 
is at the Cerenkov resonance, s = 0. Resonance at s = is possible for waves that have a component of their electric 
field along B, and this is the case for magnetoacoustic waves, but not for Alfven waves. The resonant interaction 
alone would cause the particles to diffuse in p\\ , with p^ remaining constant. This corresponds to the special case 



rf/(p) d 

dt dp\\ 



^llll(p) 



(7) 



of the general diffusion equation (A5). In the absence of pitch-angle scattering, the distribution of fast particles 
becomes increasingly anisotropic, favoring over p±, leading to a suppression of the acceleration. In the presence of 
efficient pitch-angle scattering, one averages (7) over pitch-angle, leading to the isotropic diffusion equation (3), with 
Dpp{p) = I / dfi/ji^D^w (p), where fi = p\\/p is the cosine of the pitch angle. The resulting expression for Dpp{p), given 
by (2) with (5), confirms that this resonant interaction is equivalent to Fermi acceleration. 
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III. RESONANT SCATTERING 

Effective scattering is an essential ingredient in stochastic acceleration and in diffusive shock acceleration. The 
theory of resonant scattering was developed initially in two different contexts: the diffusion of GCRs through the ISM 
and the formation and stability of the Earth's Van Allen belts. 

A. Qualitative description of resonant scattering 

Consider the gyroresonance condition (6) for s = ±1 and for waves with frequency much less than the gyrofrequency, 
<C The resonance then requires that the phase speed, cj/|/c|||, of the wave be much smaller than the parallel 
speed, \v\\\, of the particle, and the wave is equivalent to a static fluctuation in the magnetic field. The resonance 
corresponds to the wavelength, 27r/|/c|||, of the magnetic fluctuation being equal to the distance 27r|i;|||/0 that the 
particle propagates along the field line in a gyroperiod. As a consequence, the Lorentz force, qv x S'B, due to the 
magnetic fluctuation, 5B, in the wave is in the same direction at the same phase each gyroperiod. This causes a 
systematic change in the momentum of the particle, without any change in energy, corresponding to a change in the 
pitch angle, a = arccos(p||/p). The sense of the change remains the same until phase coherence is lost, resulting in 
relatively large, random changes in a. As a result, the particles diffuse in pitch angle. The resonance condition for a 
given wave (given cj, /cy) defines a resonance curve in velocity {v\\-v±) space or in momentum {p\\-p±) space, which 
is a conic section in general and which reduces to a circle in the limit u ^ 0. Pitch-angle diffusion corresponds to 
particles diffusing around the circumference of this circle. When the assumption a; ^ is relaxed, there is a change 
in the energy of the particle, with the changes in e and p\\ in the ratio u : /cy . 

The resonant interaction results in the waves either growing or damping, depending on whether energy is transferred 
from the particles to the waves, or vice versa. For any given distribution of particles, all the particles that can resonate 
with a given wave lie on the resonance curve. If the gradient in the distribution function along this curve is an increasing 
function of £, more particles lose energy to the wave than gain energy from the wave, and the wave grows. For any 
anisotropic distribution, there is a choice of the signs s = ±1, and '^y/l'^ii | that implies growth of those waves 

that carry away the excess momentum; the back reaction on the distribution of particles, called quasilinear relaxation, 
reduces the anisotropy that causes the wave growth [MW70] . This implies that the resonant waves needed to scatter 
the particles can be generated by the anisotropic particles themselves. The relevant resonant waves are Alfven and 
magnetoacoustic waves for energetic ions and for relativist ic electrons, ion cyclotron waves for nonrelativistic ions, 
and whistlers for nonrelativistic electrons. 



B. Streaming of GCRs 

Both their estimated residence time, ~ 10^ yr, in the Galactic disc and their small observed anisotropy imply that 
GCRs diffuse slowly through the Galaxy. The net flow speed is of order the Alfven speed, va ^ I0~^c. This diffusion 
is attributed to resonant scattering. A simple model for CRs streaming along the field lines at vqr is 

/(P,a) = /o(p)(l + ^cosa), fo{p) = Kck (j^^ , (8) 

where the non-streaming part of the distribution function, fo{p) is measured to be a power law above some minimum 
Po, with KcR a normalization constant. Observations imply b = 4.7 below a knee in the distribution at p ^ IQ-*^^- 
lO^^eV/c with the spectrum steepening at higher energies; Axford [A94] referred to these as GCRl and GCR2, 
respectively. 

The anisotropic CRs can generate the resonant waves that scatter them. Growth of Alfven (A) and magnetoacoustic 
(M) waves may be described in terms of a negative absorption coefficient. On evaluating the absorption coefficient 
using the general expression (A4) and inserting numerical values for cosmic rays one finds 

/ n\ o'7v.in-7/' \-'^^^fp\~^'^fcosOvcR b Va\ 

7a,m(^,6>) = -2.7 X 10 M ^) — , q— • (9) 

Vlcm -^^y \Po J v|cos^| c 3 cy 

where | cos^| is approximated by unity. It follows from (9) that the waves grow in the forward streaming direction 
(cos 6* > for vcR > 0) for -^cr > {b/3)vA' The scattering of the CRs by the waves generated by the streaming 
particles themselves causes i;cr to reduce. The rate of reduction may be evaluated using the quasilinear equation 
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written down in the Appendix. One finds 



c^'^^CR ( C - 1 |cos6>| h 

- " ^CR 7. TT o , (10) 



dt \^ 2 cosO 3 

with C = (^A ~ ^a)/(^a + -^a)' where are the fluxes of Alfven waves in the forward (+) and backward (— ) 
streaming directions, and with the scattering frequency given by 

1^3 = ^ dcosa sin^ a Daa, (11) 

with Daa determined by (A6). 

This theory works well for lower energy GCRs. The growth rate (9) is fast enough to account for growth of the 
waves that resonate with lower energy GCRs, and the scattering time is fast enough to reduce the streaming to close 
to {b/3)vA- However, the growth time increases rapidly with p, with an approximate one-to-one relation between the 
resonant wave number k and the momentum of the resonant particles kp = \q\B. For p 10^^ eV/c the growth time 
becomes ineffective, and self-generation of the resonant waves cannot account for efficient scattering. 



C. Scattering of higher energy CRs 

For higher-energy GCRs, effective scattering requires an external source of MHD waves. There is a turbulent 
spectrum in the ISM, and the observations [A95] show it to be a Kolmogorov-like spectrum 

W{k) oc (12) 

where W{k) is the energy in the waves per unit wavenumber k. Resonant scattering by a Kolmogorov spectrum 
implies a scattering frequency that decreases only slowly, Vs oc with increasing p. This scattering must 

dominate scattering by self-generated waves at higher p. It is likely that the distribution of turbulence is far from 
homogeneous, and that the most effective scattering occurs in localized regions of enhanced turbulence. However, 
how effective scattering of higher energy particles occurs is a poorly understood aspect of acceleration mechanisms. 



D. Nonresonant versus resonant growth 



Plasma instabilities can occur in two forms: resonant and nonresonant. A nonresonant instability involves an 
intrinsically growing wave, described by a complex solution of a real dispersion equation. A recent suggestion that has 
stimulated considerable interests involves an instability [BLOl, B04] that is a nonresonant counterpart of (9) [M05]. 
The potential importance of this instability is twofold. First, it can allow the lower-energy particles to generate a 
wide spectrum of waves including those needed scatter higher-energy particles. Second, its nonlinear development can 
lead to amplification of the magnetic field, which ultimately causes the instability to saturate [BLOl, B04, RKD06]. 
A similar idea has been proposed in connection with generation of magnetic fields through a Weibel instability 
associated with the relativistic shocks that generate gamma-ray bursts [WA04]. Such nonresonant growth is an 
important ingredient in the currently favored acceleration mechanism for GCRl [H06] . 



E. Formation of the Van Allen belts 



Resonant scattering was first recognized in connection with particles trapped in the Earth's radiation belts, and it 
became an important ingredient in understanding how the radiation belts are formed and maintained. Under quiescent 
conditions, ions and electrons precipitate steadily (albeit in localized bursts) from the radiation belts, implying that 
there is a continuous resupply of fast particles. The precipitation is interpreted in terms of pitch-angle scattering of 
trapped particles into the loss cone. Pitch- angle scattering of ions [W61, D61] involves ion cyclotron waves, and the 
scattering of electrons involves whistlers [D63, KP66]. 

A simple model for the formation of the radiation belts is based on adiabatic invariants, which are associated with 
quasi-periodic motions. The first adiabatic invariant is associated with the gyration about a field line, and it implies 
that M (x p\/B is an invariant. After averaging over the gyrations, a fast particle bounces back and forth along 
a field line between two reflection points. The adiabatic invariant associated with this bounce motion, denoted J, 
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FIG. 1: An idealized loss-cone distribution is illustrated; ao is the loss cone angle, the shaded region is filled by isotropic 
thermal particles, and the circular arcs denote contours of constant /. 

is the integral of along the field lines between the refiection points. After averaging over bounce motion, the 
remaining motion is a periodic drift around the Earth. There is a third adiabatic invariant, which corresponds 
to the integral for the vector potential around the Earth. For a dipole-like magnetic field it is convenient to label a 
field line by the radial distance, L, to its midpoint. Then one has M oc p^L^, J oc pyL, $ oc 1/L. Violation of the 
third adiabatic invariant, due to perturbations in the Earth magnetic field with a period equal to that of the drift 
motion, causes the particles to diffuse in L at constant M and J. The steady state density that would result from 
such diffusion is proportional to v\v\\ ^ '^/L^- Particles that enter the magnetosphere from the solar wind tend to 
diffuse inwards to build up this density. As the particles diffuse to smaller L, their momentum increases according 
to p oc L~^/^. The inward diffusion ceases at the point where pitch-angle scattering (which violates conservation of 
M) becomes effective and causes the particles to be scattered into the loss cone and to precipitate into the Earth's 
neutral atmosphere. Due to their loss-cone anisotropy, illustrated schematically in figure 1, the anisotropic particles 
can generate the resonant waves needed to scatter them, with the growth rate increasing rapidly with decreasing L. 
In this way, resonant scattering provides the sink (at smaller L) that allows the Van Allen belts to remain in a steady 
state when the magnetosphere is quiescent. The absorption coefficient for the waves tends to be maintained near zero, 
referred to as marginal stability, such that sporadic localized bursts of wave growth lead to sporadic localized bursts 
of particle precipitation. Such bursts of wave growth can be triggered by artificially generated VLF signals [H67] . 

IV. DIFFUSIVE SHOCK ACCELERATION 

Diffusive shock acceleration (DSA) is now recognized as the most important acceleration mechanism in astrophysical 
plasmas: its identification was a major achievement in the field. Prior to the identification of DSA, it was difficult to 
explain how acceleration can be as efficient as it appears to be, and why power-law energy spectra are so common. 
DSA is so efficient that one needs to look for self-regulation processes that limit its efficiency, and DSA implies 
power-law spectra of the form observed without any additional assumptions. 

A. Description of DSA 

Fermi [F54] suggested that the acceleration can be efficient if all reflections are head-on. A simple example is when 
two magnetized clouds approach each other, and a particle bounces back and forth reflecting off the approaching 
clouds. A related example is for a shock propagating into a closed magnetic loop, where a particle trapped in the 
loop can reflect head-on at the approaching shocks [W64, McL71]. However, such examples are constrained by the 
fact that acceleration ceases when the approaching structures meet. In DSA a single shock can accelerate particles 
efficiently, provided that the particles are scattered on either side of the shock [A77, K77, B78, BO], as illustrated in 
flgure 2. 

To understand this, flrst note that the fluid velocity changes discontinuously across the shock. Let the two sides of 
the shock be labeled 1 and 2. Consider a particle on side 1 about to cross the shock and enter side 2. The scattering 
centers on side 2 are convected along with the plasma, and the particle sees them approaching it head on at \ui—U2\- 
Once the particle crosses the shock and is scattered it gains energy due to the oncoming motion of the scatterers. 
After being scattered a number of times on side 2, the particle can diffuse back to the shock and cross back onto 



downstream 




FIG. 2: Diffusive shock acceleration is illustrated: the shaded vertical region is the shock, the circular blobs denote idealized 
scattering centers, and the solid line with arrows denotes the path of an idealized fast particle. The coordinate z and the 
velocities ui and U2 introduced in (14) are shown for the case of a parallel shock. 

side 1. On doing so, it sees the scattering centers on side 1 approaching it head on at \u\ —U2\. Again it gains energy 
on being scattered. DSA requires efficient scattering, and this can be achieved only by resonant scattering. Upstream 
of the shock the density of the fast particles decreases with distance from the shock, and the spatial gradient can 
cause the resonant waves to grow. Analogous to the scattering of streaming CRs, the growth of the resonant waves 
decreases rapidly with increasing particle energy, and some other source of resonant waves is required for higher energy 
particles. Scattering downstream of the shock is less problematic, with several possible sources of resonant waves, 
including waves generated in the upstream region and swept back across the shock. 

The treatment of DSA given below is a nonrelativistic, single-particle theory. The assumption that collective effects 
of the fast particles can be neglected is not necessarily valid: DSA is so efficient that the fast particles can become 
an important component of the upstream plasma. Once the pressure associated with the fast particles becomes 
comparable with the thermal pressure, the structure of the shock is modified by this pressure. This can result in the 
stresses being transferred from the downstream to the upstream plasma primarily through the fast particles, with no 
discontinuity in the density of the thermal gas [MDOl]. Such nonlinear effects provide a constraint on DSA. 

B. Diffusive treatment of DSA 

Consider a distribution of particles /(p, z) that is averaged over pitch angle and is a function of distance z from 
a shock in a frame in which the shock is at rest. It is assumed that scattering causes the particles to diffuse in the 
z direction with diffusion coefficient k{z). The particles are also assumed to be streaming with the streaming speed 
u. The diffusion is described by 



dt 



df{p, z) 
dt 

dt 



d_ 
dz 



K{z)^^^]^Q{p,z)-Uc{p), 



dz 



Jf{p.z) ^ ^df{p,z] 



dz 



dp 



P = -^P^ 



1 du 
S'd-z^ 



(13) 



where Q{p, z) is a source term, and where the sink term fesc{p) takes account of escape of particles downstream from 
the shock. The term involving a partial derivative with respect to p determines the energy changes. It is assumed 
that the speed u changes abruptly across the shock: 



ui for 2: < (upstream), 
U2 for 2: > (downstream). 



^ = {Ui -U2)S{z). 



(14) 



A stationary solution of (13) exists when both the source and the sink term are neglected, such that the equation 
reduces to udf /dz = (d / dz){ndf / dz)\ a general solution is 



/(p, 2:) = A + Bexp 



u [ dz 



(15) 
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with u constant on either side of the shock. In the upstream region, ^ < 0, one has u/k{z) ^ 0, and the solution (15) 
diverges at ^ ^ — oo unless one has 5 = for ^ < 0. Writing f±{p) = lim;2^±oo f{p^ one has 



f{p,z) 



/_(p) + [/+(p)-/_(p)]exp 
f+{p) 



Ml 



^0 



dz 



k{z) 



for z <0, 
for z > 0. 



(16) 



On matching the discontinuity in the derivative of this solution with the discontinuity due to the acceleration term, 
one finds 



Jo 



3u 



U2 - Ui 



(17) 



which determines the downstream solution in terms of the upstream solution. That is, if the shock propagates into a 
region where the distribution is /_ (p) , then after the shock has passed the distribution is /+ (p) . 

For monoenergetic injection, f-{p) oc S{p — pQ) say, (17) implies /+(p) oc p~^. The power law index, 6, is determined 
by the ratio Ui/u2, which is determined by the strength of the shock. In terms of the Mach number M = ui/vai, 
where vai is the Alfven speed upstream of the shock, and the ratio F of specific heats, one finds 



3r 3(F + 1)M^ 



r - 1 2(M2 - 1) 



(F + 1)M2 
2 + (F-l)M2' 



(18) 



where r is the compression factor across the shock. For F = 5/3 and ^ 1, the maximum strength of the shock is 
r = 4 and the minimum value of the spectral index is 6 = 4. 



C. Alternative treatment of DSA 



An alternative treatment of DSA is as follows. Consider a cycle in which a particle with momentum p crosses 
the shock from downstream to upstream and back again. Let A(p) be the change in p in a cycle, and P{p) be the 
probability per cycle of the particle escaping far downstream. Simple theory implies 

^ , , 4(r — 1) ui , 4iii 

3r V vr 

After one cycle, p ^ p' = p -\- A{p) implies that the number of particles per unit momentum changes according 
to 4:7rp^ f{p)dp A.T^p''^ f(p')dp' = [1 — P(p)]4:7rp'^ f(p)dp. After integration, one imposes the requirement that the 
downstream flux of particles balances the upstream flux, and a result equivalent to (17) is derived. 

This model allows one to estimate the acceleration time, ta, in terms of the cycle time tc- One has tc = {Xi/ui + 
^2/^2) = (Ai + rX2)/ui, where Ai^2 are the scattering mean free paths in the upstream and downstream regions, 
respectively. For isotropic scattering, the scattering mean free path. A, is related to the spatial diffusion coefficient, 
hy K = Xv/3, where v is the speed of the particle. The mean free path may be written as A = ryrg, where rg 
is the gyroradius, and it is usually assumed that ?y is a constant of order unity, called Bohm diffusion. The rate of 
momentum gain is given by 

dp_P__fdp\ + \_ 3r(Ai +rA2) 

dt't. UJioss' ^"^A' 4(r-l) ' ^'"^ 

where tg, is an acceleration time, A is a mean scattering free path, and where a loss term is included. 



D. Injection and preacceleration 

DSA requires efficient scattering. For ions, scattering by Alfven and magnetoacoustic waves requires v ^ va^ 
and for nonrelativistic electrons, scattering by whistler waves requires v ^ 43t'A- These conditions are typically not 
satisfied for thermal particles, and one needs to postulate an injection or preacceleration mechanism to create a seed 
population of ions or electrons above the respective thresholds for DSA to operate on them. 

One possible type of injection mechanism involves collisions populating a depleted thermal tail. The thresholds 
(v ^ va for ions or 1; 43i;a for electrons) are somewhere in the tail of the thermal distribution of particles, and 
the acceleration removes these particles, leaving a depleted thermal tail. Collisions between thermal particles cause 
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a net flux of particles in momentum space into the high-energy tail, tending to restore a thermal distribution. The 
rate at which nonthermal particles are created by this process depend on the collision frequency and on the details 
of the acceleration mechanism [G60]. However, this process encounters two serious difficulties. First, it is too slow 
to account for the required injection. Second, for ions it is extremely sensitive to the charge and mass, implying 
that the relative isotopic abundances of accelerated ions should differ by very large factors from the relative isotopic 
abundances of the background plasma, and this is simply not the case in general. Other preacceleration mechanisms 
involve either electric fields associated with the shock itself, or waves associated with the shock. These are discussed 
briefly below in connection with shock drift acceleration. 



E. DSA at relativistic shocks 



In the extension of DSA to highly relativistic shocks, the assumption of near isotropy of the pitch angle distribution 
is not valid, and the simple analytic theory does not generahze in a straightforward way. Numerical calculations show 
that the resulting spectrum is a power law, with b = 4.2-4.3 [HD89, AGKGOl]. 



V. DSA AT MULTIPLE SHOCKS 



In many of the applications of DSA, it is unlikely that a single shock is responsible for the acceleration [A90]. 
Extension of the theory to multiple shocks shows that DSA is related to Fermi-type acceleration: there are energy 
gains associated with DSA at each shock, and energy losses due to decompression between shocks. It is straightforward 
to include synchrotron losses in such a model, and the competition between acceleration and synchrotron losses can 
possibly explain flat or slightly inverted synchrotron spectra in some sources. 



A. Multiple shock model 



Consider the following model for DSA at multiple shocks. 

(a) All shocks are identical with compression ratio, r = 6/(6 — 3), and all have the same injection spectrum, ^0(^)5 

assumed to be monoenergetic, (/)o(p) (x S{p — po). 

(b) The distribution downstream of the first shock, fi{p), results from acceleration of the injection spectrum, (/>o(p), 

at the shock, cf. (17): 

hip) = bp-' / dqq>'-'Mq), b = -. (21) 

Jo r — 1 

(c) Between each shock a decompression occurs. Assuming that scattering keeps the particles isotropic, the distri- 

bution function after decompression is f[{p) = fi{p/R), with R = r~^l^ . Hence the injection spectrum into the 
second shock is 

rvIR 

f[{p) = b{p/R)-' / dqq'-'MQ)- (22) 

(d) The injection spectrum at any subsequent shock consists of the sum of (l)o{p) and the decompressed spectrum 

resulting from acceleration at the previous shock. 

The model implies that downstream of the nth shock (after decompression) the distribution is 

Kb' f p \ (In p/R'po)'' 



After an arbitrarily large number of shocks, the spectrum approaches /oo(p) p~^ at p > po, as illustrated in figure 3. 
This spectrum is harder, by one power of p, than the spectrum (x p~'^ from a single strong shock. Moreover, the 
spectrum /oo(p) oc p~^ is approached irrespective of the strength of the shocks, although the stronger the shocks the 
faster it is approached. An interpretation is that the combination of DSA and decompression leads to a Fermi-like 
acceleration mechanism: the asymptotic solution for Fermi-type acceleration for constant injection at p = is a 
power law with 6 = 3 for p > po- Such a distribution corresponds to a flat synchrotron spectrum: the intensity of 
synchrotron radiation is a power law in frequency, z/, with I(y) oc a = (b — 3)/2. 
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FIG. 4: The cumulative effect of DSA plus synchrotron losses after many shocks with injection at each shock; (a) the distribution, 
(b) the slope of the distribution for p > pq\ r = 3.8, pc/po = 10^, N = I, 10, 30, 50. 



B. Inclusion of synchrotron losses 



Synchrotron losses cause the momentum of the radiating electron to decrease at the rate dp/dt = —Ap^^ with 
A (X . Synchrotron losses limit DSA: there is a momentum, p = pc, above which the synchrotron loss rate 
exceeds the rate of acceleration, and acceleration to p > pc is not possible. The average acceleration rate over 
one cycle due to DSA is dp/dt = A{p)/tc = Cp, and with the inclusion of synchrotron losses this is replaced by 
dp/dt = Cp — Ap^ = Cp(l — p/pc) with pc = C / A. It is straightforward to repeat the calculation of DSA at multiple 
shocks including the synchrotron losses [MC97]. In figure 4 the logarithm of the distribution function is plotted as 
a function of log(p/po)- The synchrotron cutoff momentum is chosen to be p^ = lO^po, and all the shocks have 
r = 3.8. The distribution below the synchrotron cutoff, p < Pc, due to the cumulative effect of injection at every 
shock, becomes harder than 6 = 3, such that the slope has a peak (with b ^ 2) just below Pc- Such a distribution 
corresponds to a weakly inverted spectrum with a peak just below a sharp cutoff due to synchrotron losses. 



VI. APPLICATIONS OF DSA 



There in an extensive literature on applications of DSA, and only a few remarks on some of the applications are 
made here. The applications may be separated into astrophysical, solar and space-physics, and the emphasis in each of 
these in qualitatively different. In astrophysical applications the emphasis is on the acceleration of highly relativistic 
particles, particularly GCRs and synchrotron-emitting electrons. In solar applications one needs to explain detailed 
data on the properties of the solar energetic particles (SEPs) observed in space, and on the radio. X-ray and gamma- 
ray signatures of them in the solar atmosphere. In space-physics, the objective is to model the detailed distributions 
of particles and fields measured in situ. The comments here indicate these different emphases. 
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A. Acceleration in young supernova remnants 

Recent observational and theoretical developments have been combined in the currently favored suggestion that 
GCRs are accelerated at the shocks where young supernova remnants (SNRs) interact with the interstellar medium. 
The observational development includes X-ray, gamma ray and TeV photon data from young super novae, including 
the remnants of historical supernovae <^ 10^ yr old, that show edge-brightening, which is interpreted in terms of 
synchrotron emission from relativistic electrons accelerated at an outer shock [vA04, V06, H06]. The suggestion is 
that GCRs are also accelerated at these shocks. The theoretical development is the suggestion [BLOl, B04] that the 
CRs can cause the magnetic field to be amplified by a very large factor in association with these shocks. 

Since the original suggestion by Baade and Zwicky, it has been widely believed that the acceleration of GCRs is due 
to shocks generated by supernovae, and DSA provides an eff"ective mechanism by which this occurs. Hillas [HOG] gave 
four arguments in favor of the location being the shocks in young SNRs: (a) DSA produces a power-law spectrum of 
about the right slope, (b) The model is consistent with the long-standing arguments that the only adequate energy 
source for GCRs is supernovae, and that DSA is the only mechanism that is efficient enough in transferring this 
energy to the particles, (c) The composition of GCRl is consistent with ordinary interstellar matter being injected 
into DSA, and it is at least plausible that such composition-independent injection occurs at such shocks, (d) When 
coupled with self-generation of magnetic field, DSA can explain the energy at which the knee between GCRl and 
GCR2 occurs [H06b], cf. however [P06]. Although the details are not complete, these arguments provide a plausible 
overview as to how GCRl are accelerated. 

The amplification of the magnetic field is attributed to the nonlinear development of a nonresonant plasma instability 
driven by the CRs [BLOl, B04, RKD06]. A strong, self-generated magnetic field allows DSA to accelerate particles 
to higher energy than otherwise, due to the maximum energy being proportional to B. 

B. Acceleration in astrophysical jets 

Astrophysical jets are associated with accretion discs, around proto-stars and various classes of compact objects, 
including the super massive black holes in active galactic nuclei (AGN). Radio jets associated with AGN can extend 
to enormous distances from the AGN; these jets require that the synchrotron-emitting electrons be accelerated in situ 
[BM82]. The appearance (characterized by various knots) of the jet in M87 [MRS96] in radio and optical is remarkably 
similar, whereas the much shorter synchrotron loss time of the optically-emitting electrons leads one to expect the 
optical emission to be much more sharply localized about acceleration sites. A possible model for the acceleration 
involves DSA at multiple shocks, which develop naturally in the flow, provided the acceleration is sufficiently fast to 
overcome the losses of the optically-emitting electrons. The required scattering of the optically-emitting electrons can 
plausibly be attributed to a Kolmogorov spectrum of turbulence [M99] . 

A subclass of radio sources associated with AGB have flat or inverted spectra, corresponding to power laws oc i/^ 
with a ^ or a > over at least a decade in frequency. Synchrotron self-absorption was proposed for such spectra 
[K69] . A specific model can account for fiat spectra [M77] , but it requires such special conditions that it was referred 
to as a 'cosmic conspiracy' [C80]. A possible alternative explanation, for spectra with a < 1/3, is in terms of an 
acceleration mechanism that produces a power-law electron spectrum with b = 3 — 2a. DSA at a single nonrelativistic 
shock produces a spectrum with 6 > 4 corresponding to a < —0.5. A relativistic shock can produce a somewhat harder 
spectrum but cannot account for fiat spectra. DSA at multiple shocks tends towards b = 3, a = 0, and provides a 
possible alternative explanation for fiat spectra. The pile-up effect [S94] occurs naturally when synchorotron losses 
are included in DSA at multiple shocks [MC97], and can account for inverted spectra with a < 1/3. 

C. Solar energetic particles 

It has long been recognized that acceleration of fast particles occurs in the solar corona in connection with solar 
ffares. There is a power-law distribution of ffares in energy or in area: the energy released in a ffare is approximately 
proportional to the area that brightens in Ha. All ffares produce fast nonrelativistic (10-20 keV) electrons, observed 
through their radio (type III bursts) or X-ray (bright points) emission. Large ffares also produce energetic ions and 
relativistic electrons. In the early literature [WSW63], it was assumed that the acceleration occurs in two phases: 
the ffrst (impulsive) phase involving nonrelativistic electrons (and perhaps ions), attributed to some form of 'bulk 
energization', and the second phase involves a slower acceleration of solar CRs, attributed to the shocks that produce 
type II radio bursts. However, this simple picture was not consistent with subsequent data: gamma-ray data showed 
that the relativistic electrons and energetic ions appeared immediately after the onset of the ffare, without the predicted 



12 




FIG. 5: A cartoon [B97] illustrating one suggestion for how magnetic reconnection leads to acceleration of particles in a solar 
flare. 

delay. Spacecraft data on solar energetic particles (SEPs) revealed further phenomena and correlations that are not 
consistent with the two-phase model. Various different acceleration mechanisms have been explored for SEPs [M+97]. 

A new paradigm has emerged for the interpretation of SEPs observed in space: events are separated into flare- 
associated SEPs and CME-associated SEPs [L05]. In a CME (coronal mass ejection) a previously magnetically bound 
mass of corona plasma becomes detached, accelerates away from the Sun, and drives a shock ahead of it. The 
detachment of a CME involves some form of magnetic reconnection in the corona. Earlier ideas on a tight correlation 
between flares and CMEs have not been conflrmed by more recent data. Although there is some correlation, it is not 
one-to-one, and flares and CMEs are now regarded as distinct phenomena. The acceleration of CME-associated SEPs 
is plausibly due to DSA, but how the flare-associated SEPs are accelerated remains poorly understood. 

There are anomalies in the ionic composition in flare-associated SEPs, the most notable of which concerns ^He. 
The ratio of ^He to "^He in the solar photosphere is 5 x 10~^, but the ratio in flare-associated SEPs is highly variable, 
and can be greatly enhanced, even exceeding unity in exceptional cases. The favored suggestions for preferential 
acceleration of ^He and other isotopic anomalies is preferential preacceleration. A speciflc mechanism proposed for 
^He is the generation of ion cyclotron waves between the proton and ^He cyclotron resonances, with these damping at 
the ^He cyclotron resonance [F78]. The formulation of a self-consistent model for the acceleration of flare-associated 
SEPs is complicated by the paucity of direct signatures of processes in the solar corona, and the difficulty in reconciling 
the data on various emissions from energetic particles precipitating into the photosphere with the SEP data. 

Bulk energization of electrons is not understood. A large fraction of the energy released in a flare goes into bulk 
energization, requiring acceleration of virtually all the electrons in the flaring region of the corona . One model 
that involves shock acceleration is illustrated in flgure 5: magnetic reconnection leads to an outflow of plasma from 
the reconnection site, with the upflow associated with a CME and the downflow running into an underlying closed 
magnetic loop. In such a model, the electron heating is identifled with shock acceleration where the downflow is 
stopped by the underlying loop. There are also various non-shock-associated mechanisms, involving E\\, for bulk 
energization. 

VII. ACCELERATION BY PARALLEL ELECTRIC FIELDS 

Acceleration by an electric fleld parallel to the magnetic fleld is the simplest conceivable acceleration mechanism, 
but also the least understood. The argument for acceleration by parallel electric flelds is compelling in connection 
with aurorae, is very strong in connection with pulsars, and is somewhat controversial in connection with solar flares. 
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A. Non-MHD effects 

Acceleration by parallel electric fields, E\\, is the accepted mechanism for auroral electrons. An essential aspect of 
any acceleration by is the associated parallel current, Jy: the relevant source term from Maxwell's equations for 
the energy transfer is J\\E\\. A generic model for such acceleration requires a mechanical driver, which sets up an 
EMF, and a closed circuit around which the current flows. The energy transport away from the driver is through a 
Poynting flux, propagating at va- This energy is transferred, at least in part, to energetic particles, in an acceleration 
region, which acts as a load in the circuit. The load acts like a resistor, with the available potential localizing across 
the acceleration region, producing the . Suggested astrophysical applications of such circuit ideas tend to be viewed 
with suspicion, especially within the MHD community. This suspicion is partly justified due to some proposed circuit 
models being grossly oversimplified. However, the fact remains that there is compelling evidence for acceleration by 
E\\, and the existence of E\\ is explicitly excluded in ideal MHD. Some non-ideal- MHD effect is needed to explain 
why any E\\ develops. One such mechanism is resistivity, which is included in non-ideal or resistive MHD, but the 
E\\ allowed by classical resistivity is too weak to play any significant role in the fully ionized plasmas of interest. 
More favored possibihties include specific structures, such as double layers, or involves of Alfven waves, with the E\\ 
associated with inertial effects on the Alfven waves. 



B. Available potential 

It is sometimes helpful to note that all relevant acceleration mechanisms can be attributed to an inductive electric 
field, and to identify the specific inductive field. There are two general approaches that provide useful semi-quantitative 
estimates of the available potential. One of these is the EMF associated with changing magnetic flux through a circuit; 
the other involves the power radiated by a rotating magnet in vacuo. 

A simple expression for the EMF is in a circuit is $ = d(BA)/dt, where A is the area of the circuit. Although 
$ can be due to a changing magnetic field, $ = AdB/dt, in most cases of interest it is due to a moving circuit, 
$ = BdA/dt = BLv, where L and v need to be identified in a specific model. A potential difference appears along a 
field line that connects two regions in relative motion to each other. The magnetic field is assumed to be frozen into 
both regions (MHD applies in both), and then $ is the potential difference between them. This simple argument needs 
to be complemented with an identification of how the current closes, and how this potential is distributed around the 
closed circuit. This idea provides a plausible estimate for the energy of auroral electrons and for electrons accelerated 
in Jupiter's magnet osphere. In the case of auroral particles, the velocity, v, is between the Earth and the solar wind, 
and the acceleration occurs along the auroral field lines that connect these two regions. The potential available in this 
case can be estimated from parameters B and L across the magnetotail and v for the solar wind, giving $ of several 
kV, consistent with the observed auroral electrons of a few keV. In an early model for the lo-Jupiter interaction, the 
moon lo was assumed to drag the Jovian magnetic field lines threading it through the corotating magnet osphere, 
leading to the estimate $ = BioLioVio, with Bio is the Jovian magnetic field at the orbit of Jo, Ljo is the diameter 
of lo and vjo is the relative velocity of lo to the corotating magnetic field hues [GL69] . This gives $ of a few MV, 
which is assumed to appear around a circuit along the field lines through lo and closing across field lines in the Jovian 
ionosphere. 

A more subtle application of this idea provides a limit on the maximum energy to which particles can be acceleration 
due to DSA. The relative flow on opposite sides of the shock implies an electric field in the plane of the shock front, 
and one may interpret Bv as this electric field with v the relative flow speed perpendicular to B. Then $ = BvL is 
the potential available across the lateral extent, L, of the shock. The maximum energy to which a particle can be 
accelerated through DSA is limited to < eBvL per unit charge. 

A point of historical interest is that the earliest recognized suggestion for acceleration of CRs involved an inductive 
field due to dB/dt: in 1933 Swann [S33] proposed that the acceleration occurs in starspots due to the magnetic field 
increasing from zero to B = 0.2 T in a time r = 10^ s in a region of size L = 3x 10^ m, giving $ = BL^ /r = 2x 10^^ V. 
Although this model is not considered realistic, it is of interest that a similar value for $ results from $ = BLv with 
the same B,L and with L/r replaced by v = 3kms~-'^; these numbers are close to what one would estimate for the 
potential difference between field lines on opposite sides of a rotating sunspot. Such large potential drops are available 
in the solar atmosphere, but the implications of this are poorly understood. 

In a simple model for the electrodynamics of a rotating magnetized compact object, such as a neutron star or a 
black hole, the slowing down is estimated assuming that the rotational energy loss is due to magnetic dipole radiation. 
The power radiated may be written as where the effective displacement current, /, and the available potential, 
are related by $ = ZqI, with Zq = /xqc = 3770 the impedance of the vacuum. Equating $^/Zo to the rate of 
rotational energy loss provides as estimate of the potential available. Although the slowing down for pulsars is not 
actually due to magnetic dipole radiation, this estimate remains valid for the potential available along field lines in 
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FIG. 6: A cartoon [B96] showing the regions in the Earth's magnetosphere associated with the acceleration of auroral elec- 
trons. The generation region is where magnetic reconnection occurs, launching Alfven waves, which accelerate electrons in an 
acceleration region where the EMF from the generator is localized along field lines. 



the polar-cap region, which connect the star to the wind region outside the light cylinder. The potential may be 
interpreted as between the star and the light cyhnder, so that it leads to an along the polar-cap field lines that 
connect these regions. 



C. Acceleration of auroral electrons 



Magnetic reconnection in the Earth's magnetotail, in a so-called magnetospheric substorm, results in energetic 
electrons precipitating into the upper atmosphere, producing visible aurorae. A generic model involves a driver in the 
magnetotail, due to the solar wind sweeping field lines from the front of the magnetosphere building up the magnetic 
fiux and associated cross-tail current, relaxation through magnetic reconnection, transport of the released energy 
towards the Earth in an Alfvenic fiux at va, and closure of the cross-tail current along field lines to the ionosphere, 
across the field fines due to the Pedersen conductivity and returning to the magnetotail along nearby field fines [V02]. 
The energy fiux is partially converted into energetic electrons in an acceleration region, at a height > 10^ km, well 
above the location of the visible aurora (100 km) and well below the magnetotail, figure 6. The acceleration is due to 
an E\\ in the acceleration region, but the mechanism that sets up the E\\ is still not clear. Three ideas are discussed 
briefiy here. 

A general mechanism that produces a parallel potential difference between the ionosphere, $ion, and the magneto- 
sphere, $mag, along a given field line is associated with mirroring of thermal electrons. Let Rm be the mirror ratio, 
of the ionospheric to the magnetospheric magnetic fields. For 1 <C e^y/Te <C Rmi the so-called Knight relation 
[K73, L80, VR04] gives 

J„«i^($.„-W, ^= (2.):/4^^ . (24) 

where Te = rUeV^ is the electron temperature. The current Jy is assumed to close across magnetic field lines due to 
the Pedersen conductivity of the ionosphere, leading to a relation Jy = grad^(i;pE^), with = — grad^^ion, where 
Ep is the height-integrated Pedersen conductivity. However, the observed acceleration in substorms requires a more 
localized distribution of E^^ in the acceleration region than is implied by the Knight relation. 

In situ observations of E\\ in the acceleration region have been interpreted in a variety of ways. Localized structures 
sweeping across a spacecraft have been called electrostatic shocks [M80] and double layers (DLs) [B88, E04]. It is 
well estabfished in laboratory plasmas that DLs can develop, and models for DLs have been adapted to space and 
astrophysical plasma applications [B78, R89]. DLs are classified as weak or strong, depending on whether the potential 
difference is of order or much greater than Te/e, respectively. The formation and structure of DLs can involve two 
other ideas that have been discussed widely in the literature: anomalous resistivity and electron phase space holes 
[NOl]. However, the details of how localized DLs accelerate particles and produce the observed fiuxes of energetic 
electrons is still not clear. 

A third way in which an E^^ can arise is through the dispersive properties of Alfven waves [SOO, W02, V02]: 
the condition E\\ = is strictly valid for an Alfven wave only in ideal MHD. Inhomogeneities probably play several 
important roles. First, an inhomogeneity along the field lines can cause the Alfven speed to change, leading to refiection 
of Alfven waves, and localized trapping of their energy in an Alfven resonator [LS05, E06], between refiection at the 
top of the ionosphere and at an inhomogeneous region higher in the magnetosphere. Second, inhomogeneities across 
the field lines can lead to the Alfven waves developing large k±, such that the Alfven waves develop an E\\, becoming 
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FIG. 7: A density depletion in an auroral region of the Earth; inverted-V electrons carry a large upward current in such regions 
[Benson] . 

inertial Alfven waves [V02, W02]. This modification is significant for low- frequency Alfven waves with large k±: with 
Ae = c/cup the skin depth, the dispersion relation and the ratio of the parallel to the perpendicular electric fields 
become 

\m (l + A:iA2)V2' l^klxr ^ ^ 

The E\\ associated with inertial Alfven waves is a favorable mechanism for acceleration. The suggestion that the 
acceleration is due to the E^^ in an inertial Alfven waves is a relatively old idea [GB79], in particular for electrons 
accelerated by Alfven waves generated by lo's motion through the Jovian magnetosphere [Goertz]. This E\\, being 
a wave field, reverses sign, accelerating electrons in opposite directions, every half wave period; the wave period is 
assumed to be of order the Alfven propagation time across the Alfven resonator. The details of this acceleration are 
still not adequately understood [C02]. 

An extreme example of auroral acceleration occurs in inverted-V electrons. The upward current can essentially 
depopulate a localized auroral fiux tube of all thermal electrons, leading to a plasma cavity in which the inverted-V 
electrons are observed, figure 7. The 'inverted-V refers to the dynamic spectrum observed by a spacecraft as it crosses 
the cavity: lowest energies near the edges and maximum energy near the center. Inverted-V electrons correlate with 
auroral kilometric radiation [G74] . 

D. Acceleration in pulsars gaps 

Pulsars are rapidly rotating neutron stars with superstrong magnetic fields. There are three broad classes: ordinary, 
millisecond or recycled, and magnetars; only ordinary pulsars are considered here. These have periods P ^ 0.1 Is, 
surface magnetic fields B ^ 10^ lO^T and are slowing down at a rate P ~ 10~^^. The age, P/2P, and the surface 
magnetic field, B oc {Ppy^^, are estimated based on the model in which the rotational energy is due to magnetic 
dipole radiation. Equating the rate rotation energy is lost due to the observed slowing down to $^/Zo is thought 
to give a valid estimate for the potential available along field lines in the polar-cap region. For the Crab pulsar 
the rotational energy loss is approximately 10^^ W, implying $ 6 x 10^^ V; the available potentials in most other 
pulsars are somewhat smaller than this. In the most widely favored model, the electric field is assumed to cause a 
charge- limited fiow from the stellar surface, providing the 'Goldreich- Julian' charge density required by the corotation 
electric field. However, the Goldreich- Juhan charge density cannot be set up globally in this way, and the available 
potential develops across locahzed regions called gaps. The acceleration of 'primary' particles occurs in the gaps, and 
gamma rays emitted by these particles trigger a pair cascade that produces a secondary (electron-positron) plasma. 
This additional source of charge is assumed to screen E\\ outside the gaps. However, there is no fully consistent 
model and no consensus on the details of the structure and location of the gaps. The present author believes that 
existing models, based on a stationary distribution of plasma in a corotating frame, are violently unstable to temporal 
perturbations, resulting in large- amplitude oscillations in E\\ [L06]. 
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Although the details of acceleration by are quite different in the pulsar case from the auroral case, the underlying 
difficulties in understanding such acceleration are similar. A static model for acceleration by encounters seemingly 
insurmountable difficulties, and an oscillating model involves unfamiliar physics. 

E. Acceleration of nonrelativistic solar electrons 

A large fraction of the energy released in a solar flare goes into very hot electrons: all the electrons in a relatively 
large volume are accelerated in bulk to 10 20keV. Several models for bulk heating involving acceleration by E\\ have 
been proposed, but all have unsatisfactory features. One idea is to invoke a runaway process [H85]. Consider the 
equation of motion of an electron in the presence of an electric field and of a frictional force due to collisions with a 
collision frequency Vei^) = i^eiVe/v)^: where T4 is the thermal speed of electrons, is of the form 

dv eE fVeV 
at m V / 

It follows that an electron with speed 

are freely accelerated; Ed ''^^ called the Dreicer field. However, electron runaway sets up a charge separation that 
should quickly screen E\\ and suppress the acceleration. It also tends to cause Jy, which is opposed by inductive 
effects. Runaway on its own is not an acceptable acceleration mechanism. 

If bulk energization is due to acceleration by an E\\^ then the energies involved imply an accelerating $ of order 
10-20 kV. Two arguments suggest that much larger potential is not only available but also required. On the one 
hand, there is observational evidence that large flares tend to occur in flux loops carrying large currents, / ^ 10^^ A. 
If one identifies / with the observed current, then to produce the power of order 10^^ W released in a large flare 
requires $ of order 10^° V. With $ = 10^ V, to account for a power of 10^^ W would require a current of 10^^ A, 
which greatly exceeds the maximum possible direct current, ~ 10^^ A, that can flow in a coronal flux loop. On the 
other hand, as indicated above, Swann's [S33] old argument for a potential of order 10-*^^ V is consistent with a more 
plausible estimate for the available potential due to relative motion of magnetic foot points on the solar surface. How 
the available potential of 10^-10^^ V localizes into regions with E\\ ^ of length L ^ 10^V/£^|| is not explained 
adequately in any existing model. 

As with the auroral and pulsar cases, any model for acceleration by a static E\\ in the solar corona encounters serious 
difficulties. Although a large-scale oscillating E\\ might plausibly overcome these difficulties, there is no model for such 
oscillations, which would necessarily involve non-MHD effects. Modification of the ideas proposed for acceleration 
of auroral electrons by inertial Alfven waves in an Alfven resonator is perhaps a plausible starting point for a more 
plausible model for bulk energization. 

VIII. OTHER ACCELERATION MECHANISMS 

The three most important acceleration mechanisms are DSA, stochastic acceleration and acceleration by parallel 
electric fields. Various other mechanisms are relevant in specific applications. Examples mentioned above are the 
acceleration of energetic particles in the Earth's Van Allen belts and various preacceleration mechanisms for DSA. 
Some comments are made here on three other mechanisms: shock drift acceleration, gyroresonant acceleration and 
acceleration in current sheets. 

A. Shock drift acceleration (SDA) 

When a particle encounters a shock it is either transmitted through the shock front or reflected from it, and in 
either case it tends to gain energy [T80]. Shock drift acceleration (SDA) is attributed to the grad-B drift when the 
particle encounters the abrupt change in B: the scalar product of drift velocity and acceleration by the convective 
electric field, — gu x B is positive, implying energy gain. The energy gained by a particle depends on how far it drifts 
along the front. 



(26) 



(27) 
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FIG. 8: In the de Hoffmann-Teller frame the fluid velocities are parallel to the magnetic flelds on either side of the shock. 

In treating SDA it is convenient to make a Lorentz transformation to the de Hoffmann Teller frame. The frame 
most widely used in discussing a shock is the shock-normal frame, in which the shock front is stationary and the 
upstream flow velocity is normal to it. In this frame there is a convective electric fleld — u x B on both sides of the 
shock, and the point of intersection of B with the shock moves in the plane of the shock. In the de Hoffmann Teller 
frame this point of intersection is stationary, u and B are parallel to each other on both sides of the shock, and there 
is no electric fleld, flgure 8. (The de Hoffmann-Teller frame exists only if this velocity of intersection is less than c; in 
the opposite case, there exists a frame in which B is perpendicular to the shock normal.) In the de Hoffmann Teller 
frame the energy of particles is unchanged on crossing or being reflected from the shock. A simple theory for SDA is 
based on assuming that p\/ B \s conserved, which is approximately the case in a gyrophase- averaged sense [P81, D83], 
considering the conditions for reflection and transmission in the de Hoffmann-Teller frame, and Lorentz transforming 
to the shock- normal frame. 

Let quantities in the de Hoffmann Teller and shock-normal frames be denoted with and without primes, respectively. 
Let ^1 be the angle between ui (the shock normal) and Bi. Provided that the relative velocity between the two 
frames is nonrelativistic, the angle '0i is the same in the two frames, and the relative velocity is uq = ui tan'0i. In 
the nonrelativistic case, the pitch angles in the two frames are related by v' cosa' = vcosa + i^o, v' s'ma' = vs'ma. 
Reflection is possible only for sin^ a[ > B1/B2, and the maximum change in energy occurs for reflected particles at 
the threshold, a[ = ac, sin^ ac = B1/B2. The maximum ratio of the energy after reflection (ref) to before (inc) is for 
a particle with vi cosai = —1^0(1 — cos o^c), that is propagating away from the shock, and is overtaken by the shock. 
The maximum ratio is 



which is also the case in a relativistic treatment [K94]. For a strong shock, B2 4Bi, the ratio (28) is (2+v^)/2— = 
13.93. However, the energy ratio decreases rapidly away from the point in the momentum space of the incident particles 
where this maximum occurs [BMOl]. 

SDA accounts well for relevant observations of particles associated with shocks in the heliosphere, in particular, 
the Earth's bow shock and other planetary bow shocks. SDA produces a foreshock region ahead of any (convex) 
curved shock. The foreshock is the region ahead of the shock and behind the fleld line that is tangent to the shock, 
as illustrated in flgure 9. Fast particles reflected from the shock propagate along the fleld lines as they are swept 
towards the shock by the upstream flow. Electrons have higher speeds than ions, so that they propagate further before 
encountering the shock, producing an electron foreshock region that is much larger than the ion foreshock region. 
SDA should similarly populate a foreshock region ahead of any curved shock. 

Type II solar radio bursts are directly associated with shocks in the solar corona, and extending into the solar wind 
[NM85]. Acceleration by SDA is a favored mechanism for the acceleration of electrons that produce type II radio 
bursts [HP83, KOI]. SDA is also a potential candidate for preacceleration for DSA. Although both suggestions seem 
plausible, additional ideas are required to account for the details of the observations. 

One modiflcation of DSA, called shock surflng acceleration [USOl], involves the particles being trapped in the 
shock front for much longer than the simple theory of SDA implies. In other modiflcations, the assumption that 
the shock is a simple discontinuity is relaxed. Various additional possible acceleration mechanisms arise: due to a 
potential difference along fleld lines [JE87]; various instabilities that can lead to large- amplitude waves that can trap 
and preaccelerate particles ('surfatron' acceleration [McCOl]); and nonstationarity that involves the shock continually 
breaking up and reforming [SCD02]. 




max 



1 + (1 _ J3^/g,)l/2 

l-(l-Bi/B2)i/2' 



(28) 
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FIG. 9: The foreshock may be separated into an electron foreshock, nearer to the field line that is tangent to the bow shock, 
and an ion foreshock, nearer to the bow shock; the dashed line indicates the ill-defined boundary between these two regions, 
which actually merge continuously into each other. 

B. Resonant Acceleration 

MHD turbulence consists of a mixture of the Alfven and fast modes. Fermi- type stochastic acceleration involves 
only the fast mode, and may be interpreted in terms of Cerenkov resonance, at harmonic number s = 0. Resonances 
at s 7^ are possible and lead to dissipation of turbulence in both modes, with the energy going into the particles. For 
highly relativistic particles, resonances at high harmonics lead to a diffusion in momentum space, somewhat analogous 
to Fermi-type acceleration, but with a diffusion coefficient that depends on a higher power of p [L77] . Although this 
acceleration must occur, for say GCRs and turbulence in the interstellar medium, there is no strong case for it playing 
an important role in specific appHcations. 

C. Acceleration during magnetic reconnection 

Magnetic energy release through magnetic reconnection is known to be associated with powerful events, notably 
magnetospheric substorms and solar flares, that lead to acceleration of fast particles. There have been various 
suggestions that reconnection can lead directly to acceleration of fast particles. For example, for magnetic reconnection 
in a neutral sheet, where the magnetic field is zero on a surface and B changes sign across this surface, a test particle 
that enters the neutral sheet typically emerges with a higher energy [S65]. Although such energy gains occur, they 
do not appear to constitute an important acceleration mechanism. Both theory and observation suggest that the 
acceleration is an indirect, rather than a direct consequence of the reconnection. 

A different application of heating in current sheets has been proposed in connection with pulsar winds [C90, LKOl]. 
Beyond the light cylinder there is a wind containing wound up magnetic field lines in a nearly azimuthal direction, 
with the sign of B reversing periodically with radial distance. The wind accelerates particles at a termination shock, 
where it encounters the surrounding synchrotron nebula. An outstanding problem is that the energy fiux is Foynting 
dominated near the light cylinder, but must be kinetic-energy dominated before it reaches the termination shock. 
The suggestion is that magnetic energy dissipation in the neutral sheets separating the regions of opposite B provide 
the necessary conversion of magnetic energy into particle energy [C90, LKOl]. An alternative suggestion is that the 
number of current carriers becomes inadequate leading to the development of an accelerating field [MM96] . 

IX. FUTURE DIRECTIONS 

Three generic acceleration mechanisms are widely recognized, and there are also a number of other specific mecha- 
nisms that are important in specific applications. The pioneering works of Fermi [F49, F54] anticipated two of three 
generic mechanisms: stochastic acceleration and diffusive shock acceleration (DSA), respectively. The third generic 
mechanism, acceleration by a parallel electric field, E\\ , is the least adequately understood. 

The most important acceleration mechanism in astrophysics is DSA. It is the accepted mechanism for the accel- 
eration of GCRs and for relativistic electrons in most synchrotron sources. The favored location for the acceleration 
of GCRs is at the shocks in young supernova remnants [H06]. Despite the basic theory being well established there 
are essential details where our current understanding is incomplete and where further progress is to be expected. 
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The well-established theory is a test-particle model for DSA at a single nonrelativistic shock. Several modifications 
to the theory are known to be important: the generalizations to relativistic shocks, and to multiple shocks, and 
the dynamical effects of the accelerated particles on the shock structure. There are several aspects of the required 
efficient scattering in DSA that are inadequately understood. One is the injection problem: resonant scattering of 
nonrelativistic ions and electrons involves Alfven and whistler waves, and these require that the particles already have 
speeds v > va and v > 4:3va, respectively, and preacceleration to above these thresholds is required before DSA can 
operate. For ions, the elemental abundances of GCRs suggests an injection mechanism that is insensitive to ionic 
species, whereas for flare-associated solar energetic particles there are extreme elemental anomalies (notably ^He). 
For electrons the injection problem is more severe, and less understood. Second is the long-standing question of the 
ratio of relativistic electrons to ions [GS64]. Unfortunately, in situ data on shocks in the hehosphere provide little 
guidance on this problem: evidently the shocks in synchrotron sources have much higher Mach numbers than those 
in the heliosphere, for which there is little evidence for DSA producing relativistic electrons. A third aspect relates 
to resonant scattering of the highest energy particles. Unlike lower-energy particles, the resonant waves needed to 
scatter higher-energy particles cannot be generated by the particles themselves, and other sources of the required 
waves are speculative. A suggestion that may overcome this difficulty is that the lower-energy particles cause waves 
to grow through a nonresonant instability, which leads to amplification of the magnetic field by a very large factor 
[B04], and increases the maximum energy to which DSA can operate. Fourth, spatial diffusion across field lines is 
postulated to occur in DSA at perpendicular shocks, but is incompletely understood [AB94]. 

Despite its apparent simplicity, acceleration by E\\ is the least understood of the three generic mechanisms. An 
important qualitative point is that 7^ is inconsistent with MHD, and acceleration by necessarily involves 
concepts that are relatively unfamiliar in the context of MHD. In particular, one requires a driver that sets up an EMF, 
or available potential, and this driver is typically far removed from the acceleration region where the potential partially 
localizes to produce the accelerating . The coupling between these is best understood application is to acceleration 
of auroral electrons associated with magnetospheric substorms, but how the E\\ develops remains uncertain. Existing 
models for acceleration by E\\ in 'gaps' in pulsars and in bulk energization of electrons in solar flares have serious flaws. 
The present difficulties may be solved by assuming that the relevant E\\ is oscillating, associated with inertial Alfven 
waves in an Alfven resonator in the magnet osphere, and with large-amplitude oscillations in pulsar magnet ospheres. 
Further progress in understanding acceleration by E\^ is to be expected. 

Compared with the rapid evolution in ideas in other subfields of astrophysics, the ideas involved in understanding 
particle acceleration have developed only gradually over many decades. Further progress is certain to occur, but is 
likely to involve gradual acceptance of some of the existing ideas, rather than radical new insights. 



APPENDIX A: QUASILINEAR EQUATIONS 

The quasilinear equations are written down here using a semiclassical formalism. The waves are regarded as a 
collection of wave quanta, with energy huj and momentum ^k. Waves in an arbitrary mode, labeled M, have dispersion 
relation oj = lumO^), polarization vector eM(k) and ratio of electric to total energy, RmO^)- Their distribution is 
described by their occupation number Nm{^)- A resonant interaction occurs when the gyroresonance condition is 
satisfied: 

uj - sQ- k\\v\\ = 0, (Al) 

where s is an integer, Q = \q\B/m^ is the relativistic gyrofrequency, and /cy, v\\ are the components of k, v parallel 
to B. Resonances at s > are said to be via the normal Doppler effect, and those at s < are said to be via the 
anomalous Doppler effect. Resonances at s < are possible only for waves with refractive index greater than unity. 
The effect of a wave-particle interaction is described by the probability of spontaneous emission, wm{p,K s), which 
is the probability per unit time that the particle emit a wave quantum in the wave mode M in the elemental range 
fi^k/(27r)^. For a particle of charge q and mass m, the probability of spontaneous emission is given by 

wm{p,Ks) |eM(k) • V(k,p;s)|^(^(a;M(k) -sQ-k\\V\\), 

V(k,p;s) = {v±-i—^Js{k±R),-irjv±f^{k±R),v\\Js{k±R)), 

7 m \q\B 

An advantage of the semiclassical formalism is that the Einstein coefficients imply that the probability of stimulated 
emission and true absorption are given by multiplying this probability by wm{p,^, s)NmO^)- This allows the effect 
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on the occupation numbers of the waves and the particles to be written down using a simple bookkeeping argument. 
For the waves one finds 



d7VM(k) ^ / d7VM(k) \ 
dt [ dt 



7M(k)ArM(k), 



(A3) 



spont 



where the two terms on the right hand side describe spontaneous emission and absorption, respectively, with the 
absorption coefficient given by 
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The two forms correspond to, respectively, cyhndrical and polar coordinates in momentum space. 

The evolution of the distribution of particles due to the resonant interaction includes a term, neglected here, that 
describes the effect of spontaneous emission, and a quasilinear diffusion equation that describes the effect of the 
induced processes. In cylindrical and spherical polar coordinates, this equation is 
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respectively, with the diffusion coefficients in either set of coordinates given by 
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